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I have carried out numerical first principles calculations of the pressure depen- 
dence of the elastic moduli for several ordered structures in the Aluminum-Lithium 
system, specifically FCC Al, FCC and BCC Li, LI2 AlsLi, and an ordered FCC AlyLi 
supercell. The calculations were performed using the full potential linear augmented 
plane wave method (LAPW) to calculate the total energy as a function of strain, 
after which the data was fit to a polynomial function of the strain to determine 
the modulus. A procedure for estimating the errors in this process is also given. 
The predicted equilibrium lattice parameters are slightly smaller than found experi- 
mentally, consistent with other LDA calculations. The computed elastic moduli are 
within approximately 10% of the experimentally measured moduli, provided the cal- 
culations are carried out at the experimental lattice constant. The LDA equilibrium 
shear modulus Cn — C12 increases from 59.3 GPa in Al, to 76.0 GPa in AlyLi, to 
106.2 GPa in AlsLi. The modulus C44 increases from 38.4 GPa in Al to 46.1 GPa in 
Al7Li, then falls to 40.7 GPa in AlsLi. All of the calculated elastic moduli increase 
with pressure with the exception of BCC Li, which becomes elastically unstable at 
about 2 GPa, where Cn — C12 vanishes. 



PACS Numbers: 62.20.Dc, 81.40.Jj, 61.55.Hg 
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I. INTRODUCTION 



The ductility, light weight, and high strength of Aluminum make it the best choice for 
many areas of construction. It can, however, be improved. The addition of small amounts of 
Lithium will greatly increase the strength of the alloy with little change in the ductility [|l| . 
In particular. Young's modulus increases from 70 GPa in pure Aluminum to 94 GPa when 
the atomic concentration of Lithium reaches 32% [0]. The experimental value of Cn — C\2 
in the metastable LI2 phase of AlsLi is some 50% larger than it is in pure Al 0. 

In order to study the effect of Lithium additions on Aluminum, I used the Linear Aug- 
mented Plane Wave (LAPW) method to calculate the equation of state and the pressure 
dependence of the elastic moduli of FCC Aluminum, FCC and BCC Lithium, the LI2 phase 
of AlaLi, and an ordered FCC AlyLi superlattice. The results of these calculations may be 
used as a database for fitting approximate methods, such as the Connolly and Williams 
(CW) method the generalized Ising Hamiltonian method [^, or the embedded atom 
method . Although calculations of this type have been carried out in the Al-Li system to 
determine the phase diagram the elastic moduli were not computed at that time. 

The results of these calculations are interesting in their own right because they show 
the behavior of Al^^Li over a wide range of pressures. They also demonstrate the accuracy 
which can be achieved by the LAPW method even when the difference in energies between 
the different structures is very small (less than 1 mRy). 

The outline of this paper is as follows. In Sec. y I summarize the methods used for 
calculating elastic moduli. Sec. |T| shows how to estimate the errors in the calculations. 
The results (Cjj versus pressure) are presented in Sec. |V[ Sec. ^ contains a summary of 
the results. 



II. CALCULATION OF THE ELASTIC MODULI 



A previous paper 



HH] described the method used to determine the elastic moduli of 
intermetallic alloys from total energy calculations. I will present a brief review of the method 
here. 

The total energy calculations are carried out using a full potential version [|ll|] of the 
Linear Augmented Plane Wave Method [Q. I will refer to this as "the LAPW method". 
The transformation from the many-body problem to a single-particle picture is obtained via 
the Hedin-Lundqvist parameterization |jl2| of the Local Density Approximation [T^ to the 
density functional theory |T^. The LAPW method treats the core states (here the Al Is, 
2s, and 2p states, as well as the Li Is states) fully relativistically, and treats the conduction 
bands in the semi-relativistic approximation |]T5[. The lattice is divided into muffin-tin 
spheres and an interstitial region. Inside the muffin-tins the basis functions are expanded 
into spherical harmonics up to order / = 8. In the interstitial region the basis functions are 
expanded into plane waves. I experimented with the basis set cutoff, and found that at the 
equilibrium volume approximately 40 wave functions/atom were sufficient in the sense that 
the predicted elastic moduli did not change when I increased the number of basis functions. 
The single particle potential is also also expanded into spherical harmonics and plane waves, 
however the spherical harmonic expansion inside the muffin-tins is only carried out up to 
terms of order / = 4. The Brillouin-zone integrations were performed using the Monkhorst 
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and Pack special k-point prescription [|T^], with modifications to properly treat the special 
k-points in reduced symmetry lattices [|T0[] . 

The elastic moduli of a cubic crystal may be divided into two classes, the bulk modulus 
B = (Cii + 2Ci2)/3, and the two shear moduh, Cu — C12 and C44. The bulk modulus is 
related to the curvature of E{V), 

B{V) = -VP'{V) = VE"{V) , (1) 

where V is the volume of the unit cell, E{V) is the energy/unit cell at volume V, and 
P(y) = —E'iy) is the pressure required to keep the cell at volume V. Since the calculations 
only provide a set of energies E{Vi) for a limited number of volumes Vi, the second derivative 
E"(y) must be approximated. Here I begin by making a least squares fit of the computed 
energies to the form proposed by Birch 



Q T/ ^ T/ 

- miiyr^' - 1? + e 7n[(^)^/^ - 1]- . (2) 

where Eo, Vo, Bo, and B'^ are, respectively, the equilibrium energy, volume, bulk modulus, 
and pressure derivative of the bulk modulus, while N is the order of the fit. For a second 
order fit (A^ = 2) it is obvious that B'^ = 4. Experimentally, B'^ is usually between 3 and 5. 
The bulk modulus can then be obtained by analytic differentiation of (^. 

The shear moduli require knowledge of the derivative of the energy as a function of a 



lattice strain |]T8[. In the case of a cubic lattice, it is possible to choose this strain so that 
the volume of the unit cell is preserved. The strain can also be chosen so that the energy 
is an even function of the strain, whence an expansion of the energy in powers of the strain 
contains no odd powers. Thus for the calculation of the modulus Cu — Cu I used the 
volume-conserving orthorhombic strain tensor. 



'6 \ 



(3) 



e= \ -6 

5^(1 -5^) J 

Application of this strain changes the total energy from its unstrained value to 

E{S) = E{-6) = E{0) + {Cu - Cu)V6^ + 0[6^] , (4) 

where V is the volume of the unit cell and E{0) is the energy of the unstrained lattice 
at volume V. For the elastic modulus C44, I used the volume-conserving monoclinic strain 
tensor 

¥ \ 

00. (5) 

(57(4 - 52) ) 

which changes the total energy to 

E{6) = E{-6) = E{0) + ^CuV6^ + 0[6^] . (6) 
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191 in either or (|]) since the strains @ 



Note that there is no pressure or stress term 
and d^) are constructed so that AV = 0. 

The strains (|^) and (|) can be used for any cubic lattice. In the general case, the 
internal parameters of the lattice must be chosen to minimize the total energy of the strained 
structure. Fortunately, the lattices discussed here fall into a restricted subset of the cubic 
lattices, where all of the atoms sit at inversion sites, even under the reduced symmetry 
caused by the strains (^) and (^). For lattices in this class, the atoms will remain on the 
inversion sites for infinitesimal strains 5. The sites may become unstable at some finite 
strain, but even then the force on the atoms at the inversion sites will be zero. Since the 
elastic modulus is only concerned with the limit 5 ^ 0, it is unnecessary to relax the internal 
parameters in the strained lattice. This enormously simplifies the calculation. Furthermore, 
the elimination of the linear term ensures that S = will be a extremal point, as required 
by symmetry. 

The use of the strains (^) and (^) reduces the symmetry of the problem compared to 
the tetragonal and trigonal distortions which are sometimes used ||10[- The lower symmetry 
means that more k-points are generated, however this is compensated by the fact that 
and (^) are even functions of the strain 6, and so we need only one-half as many calculations. 

Having outlined the methods used to calculate the Cjj, it must be noted that many 
compounds, including some of those in the Al-Li system, cannot be grown in single crystal 
form, hence the individual elastic moduli cannot be measured experimentally. In these cases 
experiments can only determine the isotropic bulk B and shear G moduli of polycrystalline 
aggregates of small crystallites. In the case of the cubic systems discussed here the calculated 
Gij can be used to determine B exactly and to place rather strict bounds on G. General 
bounds on B and G were originally determined by Reuss [20| and Voight |21]. Later, Hashin 
and Shtrikman found improved bounds specific to cubic materials 



We will use the later 

bounds here. For isotropic polycrystalline aggregates of cubic crystallites, the bulk modulus 
is given exactly by 
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B = -{Cn + 2Cu) 



(7) 



just as for a cubic crystal p3|. In true isotropic materials, the shear modulus is related to 
the elastic moduli by 



G' 



'-^44 



(8) 



but in real crystals the anisotropy constant 



A 



2G 



44 



(9) 



Cn — Ci2 

is not unity. In this case, we can only bound the shear modulus of the aggregate. Hashin 



and Shtrikman ||22| found that for cubic crystals these bounds are given by 

HGl-Gl) 



Gi — G* + 



5-A(3r{G*,-Gl) 



(10) 



and 
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where 



and 



- + 5 - QHG\ - Gl) ' ^^^^ 



^1 ~ 2^^^^ ~ ^12)) — C44 (12) 



_ 3(i? + 2Gt) _ 3(i? + 2G;) 

56-1(35 + 46^) ' 5G^(35 + 46^) ' ^ ' 

The Shtrikman bound Gs is designated as the smaller of Gi and G2, while the Hashin bound 
Gh is the larger. Note that in the limit of an isotropic lattice (A = 1, or = G2) we get 
Gs = Gh = G. 

In the Al-Li system these bounds are so tight that the difference \Gh — Gs\ is on the 
order of the uncertainty in the calculation of the individual Gij. Thus I will use 

G = ]^{Gs + Gh) (14) 



as the definition of the shear modulus. The estimated error in G is on the order of the 
estimated error ii 
Young's modulus 



estimated error in the individual Gij. Associated with the bulk and shear modulus are 



and Poisson's ratio 



(16) 

for isotropic crystals. Since these quantities are often of interest I will list them in the 
results. 



III. ESTIMATING ERRORS 

Using equations (|) and (|^), the elastic moduli Gu — Gu and C44 are simply related to 
the second derivative of E{6) at zero strain. If S' is an infinitesimal strain, then, since we 
have eliminated terms linear in 6, 

a = hm ^ ^ ^ , 17) 

where a is the linear combination of the Gij associated with the strain 6. In principle, we 



could pick a small, finite strain 5, calculate E(6) and E{0), and use (Jl7\} to estimate the 
elastic modulus. In practice, however, numerical errors, which cause fluctuations in the 
computed E{6), make the calculated value of a change rapidly and unpredictably with S 
unless great care is taken to ensure convergence in both k-points and basis set size. 
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To avoid this difficulty, it is better to chose a set of M points 5i and fit the resulting E{5) 
to a polynomial of the form 

N 

eiv[{/4;'^] = E/n'^"^- (18) 

n=0 

The polynomial (|1^) contains no odd powers in 5 because these terms have been eliminated 
from the expansions @j and (1). The quantity fi/V plays the role of a in Eq. (|T7D . Note that 
/i may still be very sensitive to the choice of N and the We need a reliable procedure for 
determining the best value of /i from the available data. This procedure should determine 
the order 2N of the polynomial expansion ([l8|), and also give some estimate of the possible 
values of /i which would give a nearly equivalent fit. To do this analysis, I have chosen to 



follow a method outlined in the Numerical Recipes book . 

The calculation of an elastic modulus Cij from total energy calculations begins by choos- 
ing a set of M points 5j, {i = 1,2, ...,M}. Typically I choose M = 5, with Si = and 
^5 ~ 0.05. If the 6i are equally spaced, then for the Aluminum-bearing compounds discussed 
here AEi = E{6i+i) — E{6i) is on the order of 1 mRy. Smaller values of AE may cause 
excess numerical noise to introduce large uncertainties in the analysis. 

To fit ([T8|) there must be an energy Ei and an error estimate cTj for each Si. I estimated 
Ei and cTj by performing LAPW total energy calculations on several k-point meshes. For 
FCC Li and Al I used meshes corresponding to 32000, 42592, and 55296 points in the full 
Brillouin zone. Symmetry reduced the actual number of points to no more than 7200 in the 
irreducible part of the Brillouin zone. For structures with larger unit cells and hence smaller 
Brillouin zones, correspondingly smaller numbers of k-points were used. In retrospect, this 
large number of k-points may be excessive, but AE{6) in Li is rather small (the energy 
change between 6 = and 6 = 0.1 is about 1 mRy), so high accuracy was needed. In 
any case, since the secular equation was small (about 40 basis functions), these calculations 
could be done quickly even for large numbers of k-points. The energy Ei was calculated by 
taking an average of the computed E{6i) over the different meshes, weighted by the number 
of points in the mesh. The uncertainty in the calculation, cTj, was assumed to be the k-point 
weighted standard deviation about this average. This method of determining cxj does not 
include errors associated with other variables in the calculation, including basis set size and 



the Fast Fourier Transform mesh size |2^. However, it should give at least an order of 
magnitude estimate of the size of the error in the calculation. 

For a given N, the parameters /„ in ([TsD are chosen by a least squares fitting procedure, 
minimizing the quantity 



x'{N) = f: ( EizJi^MhM] ' . (19) 



1=1 



If the numerical errors are normally distributed, the probability that an error less than or 
equal to x^(^) can occur by chance is 

q(N) = / tiM-N-3)/2^-t ^f/YUM -N- l)/2) . (20) 

For a given compound and lattice constant, the best N is the one which maximizes q{N). 
Typically, the best q was on the order of 0.5, indicating that our error estimates cTj were prob- 



6 



ably too large p4[. (Usual values of o are on the order of 0.03 mRy for these calculations.) 
Maximum values of q ranged from 0.05 to 0.95. 

As an example of this process, Fig. |1] plots the energy of LI2 AlsLi as a function of 
the orthorhombic strain for the fixed volume associated with the cubic lattice constant 
a = 7.50 a.u. (r^ = 2.93 a.u.). The = 1 and N = 2 fits are also shown. Obviously the 
A^ = 1 fit (g = 0.002) is worse than the N = 2 {q = 0.99) fit. On the other hand, as seen in 
Fig. ^, A^ = 1 (g = 0.68) and N = 2 [q = 0.74) fit the monoclinic strain energy (|^) equally 
well. 

Once the order N of the fit (|T^) has been determined, the accuracy of f\, the parameter 
related to the elastic modulus, can be established. Assume that the numerical errors in the 
calculation are normally distributed. Then there is a 68.3% probability that the correct /i is 
within the set of all possible values of /i which change x^(^) by less than 2.3. The bounds 
I present for the elastic moduli are those selected by this formula. 

The equilibrium parameters and confidence limits derived from the Birch fit {Eg, Vo, Bg, 
and -Bq) are obtained by similar methods using the functional form (0) and first-principles 
data E{Vi). 

1 have written out the above procedure to show how an estimate of the errors in these 
calculations can be obtained and so that other researchers will be able to obtain the same 
elastic moduli and estimated uncertainties from the same set of data. The accuracy of 
these uncertainties is obviously open to question. The errors here are probably not normally 
distributed. In fact, some systematic errors, such as the change in the energy with respect 
to basis set size, are not included in the analysis. To conclude, then, these error estimates 
only serve as a guide to the uncertainty in these calculations. 



IV. RESULTS 



The results are presented as a series of tables and figures. I first used the LAPW method 
to calculate the total energy versus volume for FCC Al, FCC Li, LI2 AlsLi, and ordered 
Al7Li in an FCC supercell, as diagrammed in Fig. |^. For completeness, calculations were 
also performed for BCC Li. The data from these calculations is shown in Fig. ^, where 
the energy/atom is plotted versus the Wigner-Seitz radius. The data was fit to the Birch 
equation (Q), as described in Sec. |ITJ In all cases the A^ = 3 Birch fits had the best "q" value 
(^). Table H lists the calculated equilibrium lattice constants, bulk moduli, and pressure 
derivatives for the A^ = 3 Birch fits. I also list a theoretical binding energy, calculated by 
finding the difference between the LAPW total/energy atom and the self-consistent total 
energy of a spherically averaged atom using the Hedin-Lundqvist approximation to the 
LDA and the semi-relativistic approximation [0 for the "valence" orbitals. The computed 
results are in good agreement with the calculations of Guo et al. [26,^|, which I also present 
here. 



The available experimental data for the lattice constants ||28| , ^ and bulk modulus ||30|j3T 



are also listed in Table |. The predicted lattice constants are all smaller than experiment. 
The error is largest in BCC Lithium, where the lattice constant is 3% smaller than experi- 
ment. The Al and AlsLi lattice constants are within about 1% of experiment. While some 
of this is due to the neglect of zero point motion and thermal expansion, much of the error, 
especially in Lithium, can be attributed to problems with the LDA 
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The equation of state, P{V), for these curves can be determined by differentiating the 
fit (H). The volume dependence of the bulk modulus, B(y), can be found by applying the 
formula (|1|) to Since P{V) is a monotonic function over the range of volumes shown 
here it is trivial to numerically invert this function and thus calculate B{P), the pressure 
dependence of the bulk modulus. Fig. |^ plots B(P) for the systems studied in the paper. 

The bulk moduli are in reasonably good agreement with experiment. The predicted 
Lithium modulus is about 26% too large, while the Al modulus is 8% too large. The latter 



error is typical of other calculations |]T0|. Since the LDA underestimates the equilibrium 
lattice constant, it is not surprising that it overestimates the bulk modulus. If we use the 
Birch fit to evaluate the bulk modulus (|1]) at the experimental lattice constants, we 
get bulk moduli of 10.6 GPa for BCC Lithium and 66.7 GPa for Aluminum. These values 
are both 12% lower than the respective experimental values, overcorrecting the experimental 
values. It is not clear whether this remaining error is due solely to the LDA or if contributions 
to the bulk modulus from the thermal and zero-point motion should be included. 

The shear moduli were somewhat more difficult to obtain, since the small variation in the 
total energy as a function of the shear 6 required a large k-point mesh. The uncertainties 
in these moduli are usually larger than the uncertainties in Bq. Fig. ^ and Fig. |^ show 
Cu —Ci2 and C44, respectively, as a function of pressure. The error bars show the estimated 
uncertainty in the moduli (see Sec. PTT]) . The moduli were actually calculated at fixed volume. 
The pressure was calculated from the Birch fit @) using PiV) = —E'{V). 

As one can see from the figures, the pressure dependence of the elastic moduli can be fit 
by a straight line, at least over the range —2 GPa to 12 GPa. The lines shown in the figures 
are determined by a least squares fit to the computed elastic moduli, using the estimated 
uncertainty as a weighting factor. The intercepts of these lines with the P = axis are listed 
in Table 0, while the slopes are in Table The estimated errors in the fit are calculated 
in the manner of Sec. |T|. 

Table ^ also lists the available experimental data for Aluminum |^, and the AlsLi 
calculations of Guo et al. |Q. There is a slight volume discrepancy between the computed 



Aluminum moduli and the experimental numbers, so it is not surprising that the calculated 
moduli are larger than experiment. Evaluating the moduli at the experimental volume gives 
Cii — C12 = 48.7 GPa and C44 = 30.1 GPa, both within 8% of experiment. Combining these 
numbers with the bulk modulus data shown above, it seems that the LDA can predict the 
elastic moduli to within about 10%, if the experimental volume is used. 

Except for Cu, the computed AlsLi moduli are somewhat smaller than the calculations 
in Ref. ||3^ or the experimental results reported there. Those calculated moduh were pre- 
sumably determined at the experimental volume, so the true discrepancy is somewhat larger 
than the 42% disagreement in C44 shown in the table. However, both sets of calculations 
agree that there is a large increase in Cu — Cu with increasing lithium content. 

I next used the zero pressure moduli from Table |1| to calculate the Hashin-Shtrikman 
bounds on the shear modulus G for a polycrystalline aggregate. Notice that, except for Li, 
the bounds on G are smaller than the uncertainties in the original elastic moduli. Thus 
the actual uncertainties in G are approximately the same as the uncertainties in Table ||. 
Table ^ collects the information which can be derived for the aggregate crystals, including 
the bulk modulus (|^, the shear modulus (0), Young's modulus (p!5D, and Poisson's ratio 
The Young's modulus for several Al-Li compounds has been measured experimentally 
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0, and is also presented in Table ||. Not surprisingly, the predicted Young's modulus is 
some 23% larger than the experimental number. Evaluating the modulus at the experimental 
parameter gives E = 103 GPa, within 13% of experiment. 

Although not part of the major thrust of this paper, there are several features of interest 
in the Lithium data. Calculations show that FCC Li is 0.2 mRy more stable than BCC Li. 



Experimentally, BCC Lithium may transform to the FCC state at T = 78K ||28[, although 
the actual ground state of Lithium is apparently a hexagonal close packed structure. The 
calculations are also in agreement with earlier works |0, which agree that the FCC-BCC 
energy difference is about 0.2 mRy. 

The behavior of Lithium under pressure is also worthy of note. Fig. ^ shows that BCC 
Lithium becomes elastically unstable at about 2 GPa, when the modulus Cu — Cu vanishes. 



This is consistent with the ground state of Lithium being a close-packed phase [^]. However, 
the Birch fits to the equation of state show a pressure induced phase transition from FCC 
Lithium to BCC Lithium near 0.6 GPa. Since the energy difference between the two phases 
is only 0.2 mRy, it is obvious that thermal and zero point motion, which are neglected here, 
must be included to obtain an understanding of the behavior of Lithium. 



V. SUMMARY 

As shown above, addition of Lithium into FCC Aluminum increases both of the shear 
moduli. Cu — Cu is 28% larger in AlyLi than in pure Al, and 79% larger in AlsLi than 
in Al. For C44 the increases are 20% and 6%, respectively. Conversely, the bulk modulus 
decreases slightly with increasing Lithium concentration. Since Cu — C12 and C44 increase 
by nearly the same ratio in going from Al to AlyLi, the anisotropy factor A is nearly 
unchanged, while it is much smaller in AlsLi than it is in either Al or AlyLi. This increase 



in anisotropy was also found by Guo et al. |33|, who attributed it to increasing anisotropy 



of the chemical bonding with the addition of Lithium. From the results presented here, 
it is apparent that the anisotropy ratio does not change appreciably for small amounts of 
Lithium, increasing only when there are large numbers of Lithium atoms on next-nearest 
neighbor sites (see Fig. However, Poisson's ratio ( [I6| ) does decrease from 0.32 in Al to 0.26 
in Al7Li, indicating that some anisotropy is present even at this Li concentration. It would 
be interesting to study the behavior of the anisotropy at smaller Lithium concentrations by 
continuing the elastic moduli calculations for ordered supercells with chemical composition 
AlisLi or Al26Li. Unfortunately, while simple total energy calculations for these structures 
are possible on present computers, in this case the not all of the Aluminum atoms are 
at inversion sites, so relaxation of these atoms around the Lithium "impurities" must be 
included p5|]. It would require a large number of total energy calculations to obtain the 



energies of the relaxed structures needed for the calculation of the Cij in these systems. 

The elastic moduli all change linearly with pressure in the region studied (-2 to 12 GPa). 
There is some variation in the slope of these lines, but generally the slope is between 3 and 
5. The exceptions are C[i — C'12 and C44 in the Lithium structures, which are all less than 
1, and C[i in the Aluminum compounds, where the slope is about 7. 

In conclusion, I have calculated the pressure dependence of the elastic moduli for several 
ordered Al-Li lattices, including error estimates in the calculation of the elastic moduli. The 
computations show that the addition of Lithium strongly increases the shear moduli and 
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Young's moduli (by about 20%) over Aluminum. The decrease in Poisson's ratio indicates 
that the Lithium introduces anisotropic chemical bonding. 
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FIGURES 



FIG. 1. Energy as a function of the square of the orthorhombic strain @, used to determine 
Cii — Ci2 in LI2 AlsLi. The unstrained cubic lattice constant is a = 7.50 a.u. (r^ = 2.93 a.u.). 
The error bars, representing the estimated error cjj, defined in the text, are not shown as they are 
smaher than the heights of the symbols. The solid line is the = 1 fit ([T^), the dashed line is the 
iV = 2 fit. 



FIG. 2. Energy as a function of the square of the monoclinic strain (|5|), used to determine C44 
in LI2 AlsLi. The unstrained cubic lattice constant is a = 7.50 a.u. (r^ = 2.93 a.u.). The error 
bars on the data points, which represent the errors o"j defined in the text, are omitted as they are 
almost the same height as the symbols. The solid line is the = 1 fit ([T^), the dashed line is the 
A^ = 2 fit. 



FIG. 3. The FCC based structures used in this paper. The FCC lattice sites are represented 
by open, shaded, and filled circles. When all sites contain the same type of atom, this represents 
a simple FCC lattice, with the cube volume equal to four times the volume of the unit cell. If the 
filled sites are occupied by Al atoms and the open and shaded sites by Li atoms, this represents 
AlsLi in the LI2 structure, with the cube shown being the entire unit cell. Finally, if the shaded 
and filled sites are occupied by Al atoms and the open sites by Li atoms, this represents a part of 
the FCC supercell Al/Li lattice. The eight cubes which touch at the lower left-hand corner of the 
lattice form a cube enclosing four times the volume of the AlyLi unit cell. 



FIG. 4. Computed energy/atom (in Rydbergs) versus the Wigner-Seitz atomic radius (in 
atomic units) for compounds in the Al-Li system. The energy is relative to the atomic energies 
calculated for spherically averaged atoms. (See the text.) The corresponding A^ = 3 Birch fits @ 
are plotted as continuous lines. In the upper part of the graph, the solid line and the "+" represent 
BCC Li, the dashed line and the "o" FCC Li. In the lower part of the graph, from the bottom 
up the curves are FCC Al (solid line and "□"), FCC AI7U (dashed line and "x"), and LI2 AlaLi 
(dotted line and "A"). The error bars representing the would be smaller than the size of the 
markers and so are not shown. 



FIG. 5. Computed bulk modulus versus pressure. All units are GPa. Since this is computed 
from (2), using the Birch fit (|2|) in Fig. no data points are shown. In the lower part of the graph 
the solid line represents BCC Li, the dashed line FCC Li. In the upper part of the graph the solid, 
dashed and dotted lines represent FCC Al, FCC AlyLi, and LI2 AlsLi, respectively. Over this 
pressure range B{P) can be reasonably approximated by a straight line. 



FIG. 6. Shear modulus Cn — C12 versus pressure. All units are GPa. In the lower part of the 
graph, the solid line and the "+" represent BCC Li, the dashed line and the "o" FCC Li. In the 
upper part of the graph, from the bottom up the curves are FCC Al (solid line and "□"), FCC 
AlyLi (dashed line and "x"), and LI2 AlsLi (dotted line and "A"). The calculation of the error 
bars is described in Sec. III. The straight lines are least squares fits to the calculated pressure 
dependence of Cn — C12, weighted by the estimated error. 
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FIG. 7. Shear modulus C44 versus pressure. All units are in GPa. In the lower part of the 
graph, the solid line and the "+" represent BCC Li, the dashed line and the "o" FCC Li. In the 
upper part of the graph, from the bottom up the curves are FCC Al (solid line and "□"), LI2 AlsLi 
(dotted line and "A"), and FCC AlyLi (dashed line and "x"). The calculation of the error bars 
is described in Sec. |T|. The straight lines are least squares fits to the calculated C44(P), weighted 
by the estimated error. Note that the Al7Li and AlsLi curves are reversed compared to Fig. ^. 



14 



TABLES 



TABLE L Equilibrium constants for cubic structures in the Al-Li system, using the = 3 
Birch fit (Q). The Wigner-Seitz radius is in atomic units, the bulk modulus Bq in GPa, and the 
cohesive energies Ec in Rydbergs. 

-0.1511 
-0.1245 

-0.1513 
-0.1250 

-0.308 
-0.295 

-0.292 
-0.283 
-0.277 
-0.272 





System 


rs 


Bo 




BCC 


Li^ 


3.13 


15.1±.2 


3.1±0.5 


BCC 


Li^ 


3.13 


15.0 




BCC 


Li 


3.24^= 


12.0^^ 




FCC 


Li'" 


3.16 


15.0±.l 


3.4±0.2 


FCC 


Li*' 


3.12 


14.0 




FCC 


Li 


3.24= 






FCC 


AP 


2.95 


82.4±.9 


4.8±0.1 


FCC 


AF 


2.95 


82.0 




FCC 


Al 


2.99= 


76.3^ 




FCC 


AlrLi^ 


2.94 


75.6±.6 


5.0±0.3 


FCC 


AlrLi^ 


2.94 


74.0 




LI2 


AlgLi^ 


2.93 


69.1±.5 


4.7±1.2 


LI2 


AlsLi^ 


2.94 


72.0 




LI2 


AlsLi 


2.95s 







^Calculated, this work. 
^Calculated, Ref. 



^Experiment, Ref. |28] 
'^Experiment, Ref. |3C] 
^Calculated, Ref. 
^Experiment, Ref. pTf, 



^Experiment, Ref. |29| 



TABLE II. Equilibrium elastic moduli for cubic structures in the Al-Li system. The 
Wigner-Seitz radius is in atomic units, and represents the equilibrium lattice constant for com- 
putations and the room-temperature lattice constant for experimental measurements. All moduli 
are in GPa. The last column lists the anisotropy ratio, A (M). 





System 


rs 


B 


Cii — C12 


Cii 


C12 


C44 


A 


BCC 


Li^ 


3.13 


15.1±.2 


0.8±0.1 


15.6±0.3 


14.8±0.3 


11.1±0.2 


27.8 


FCC 


Li^ 


3.16 


15.0±.l 


3.9±0.6 


17.6±0.5 


13.6±0.3 


9.1±0.2 


4.6 


FCC 


Al^ 


2.95 


82.4±.9 


59.3±1.0 


121.9±1.6 


62.7±1.3 


38.4±3.0 


1.3 


FCC 


Al^ 


2.99 


76.3 


46.0 


107.0 


61.0 


28.0 


1.2 


FCC 


AlrLi'^ 


2.94 


75.6±.6 


76.0±3.0 


126.2±2.6 


50.3±1.6 


46.1±2.7 


1.2 


LI2 


AlgLi'^ 


2.93 


69.1±.5 


106.2±3.8 


139.8±3.0 


33.7±1.7 


40.7±2.2 


0.8 


LI2 


AlgLi^ 




72.3 


128.6 


158.0 


29.4 


57.7 


0.9 



^ This work. 

^ Experiment, Ref. pl| . 

= Calculations, Ref. |33] 
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TABLE III. Equilibrium pressure derivatives of the elastic moduli for cubic structures in the 



Al-Li system. The Wigner-Seitz radius is in atomic units. 





System 




B' 


Cn C'12 


C'u 


C12 


C44 


BCC 


Li 


3.13 


3.1±0.5 


-0.3±0.03 


2.9±0.6 


3.2±0.5 


0.7±0.1 


FCC 


Li 


3.16 


3.4±0.2 


0.4±0.4 


3.6±0.5 


3.2±0.3 


0.5±0.1 


FCC 


Al 


2.99 


4.8±0.1 


3.3±0.3 


7.0±0.4 


3.7±0.2 


2.6±1.0 


FCC 


AlrLi 


2.94 


5.0±0.3 


2.2±0.7 


6.5±0.7 


4.3±0.5 


3.1±0.6 


LI2 


AlgLi 


2.93 


4.7±0.2 


4.4±1.2 


7.6±0.9 


3.2±0.5 


2.5±0.7 



TABLE IV. Derived properties of aggregate crystals. The moduli B, Gs, Gh, G, and E are in 
GPa. All calculations are done using the LDA equilibrium moduli in Table 





System 


B 


Gs 


Gh 


G 


E 


a 


BCC 


Li 


15.1 


1.6 


5.1 


3.3 


9.3 


0.40 


FCC 


Li 


15.0 


4.6 


5.4 


5.0 


13.5 


0.35 


FCC 


Al 


82.4 


34.6 


34.6 


34.6 


91.2 


0.32 


FCC 


Al7Li 


75.6 


42.7 


42.7 


42.7 


107.7 


0.26 


LI2 


AlsLi 


69.1 


45.3 


45.3 


45.3 


111.5 


0.23 


LI2 


AlgLi^ 










91.0 





Experiment, Ref. interpolated to 25% Lithium content. 
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